Transport Properties of the Lorentz Gas 
in Terms of Periodic Orbits 



PREDRAG CVITANOVIC 

Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

JEAN-PIERRE ECKMANN 

Departement de Physique Theorique, Universite de Geneve, CH-121 1 Geneva 4, Switzerland 

and 



^ PIERRE GASPARD 



Service de Chimie-Physique, Universite Libre de Bruxelles, Campus Plaine, B-1050 Brussels, Belgium 

(Received 17 October 1993) 



o 

Q 

in 

Abstract — We establish a formula relating global diffusion in a space periodic dynamical system to cycles in the elementary 
cell which tiles the space under translations. 
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The Effusive properties of a "Lorentz gas" [L] have been studied extensively in the literature, see refs. [S, BS, 
Bu, ON, MZ]. The novelty of the approach presented here is that it provides an explicit connection between 
the global diffusion and the dynamics restricted to an elementary cell. Our method applies to any hyperbolic 
dynamical system that is a periodic tiling M = \J heT M h of the dynamical phase space M by translates M h of 
an elementary cell M, with T the Abelian group of lattice translations. Furthermore, each elementary cell may 
be built from a fundamental domain M by the action of a discrete (not necessarily Abelian) group G. 

These concepts are best illustrated by a specific example, a Lorentz gas based on a Sinai billiard [S] on the 
hexagonal lattice as in Fig. 1, with disks sufficiently large so that no infinite length free flight is possible. 

It should be stressed that M refers here to the full phase space, i.e., both the spatial coordinates and the 
momenta. The spatial component of M is the complement of the disks in the whole space. We shall relate the 
dynamics in M to diffusive properties of the Lorentz gas in M, using functional determinants and ^-functions. 

It is convenient to define a time evolution operator for each of the 3 cases of Fig. 1. x t = <p (£) denotes 
the point in the global space M obtained by the flow in time t. x t = </>*(#) denotes the corresponding flow in 
the elementary cell; the two are related by 



n t (x) = 4>\x) - ^(x) G T , 



(1.1) 
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Fig. 1. Tiling of M, a periodic lattice of reflecting disks, by the fundamental domain M. Indicated is an example of a global trajectory 
x t together with the corresponding elementary cell trajectory x t and the fundamental domain trajectory x t . 

the translation of the endpoint of the global path into the elementary cell M. The quantity x t = denotes 
the flow in the fundamental domain M; </>*(£) is related to </>*(£) by a discrete symmetry g £ G which maps 
x t £ M to x t G M. 

Fix a vector (3 G R d , where (Z is the dimension of the phase space. We will compute the diffusive properties 
of the Lorentz gas from the generating function 



where the average is over all a; £ M. 

The diffusive properties follow by studying 
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and its derivatives at /? = 0. Clearly Q(0) = 0, and if by symmetry all odd derivatives vanish, there is no drift 
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In that case, the second derivatives 
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yield a (generally anisotropic) diffusion matrix. The spatial diffusion constant is then given by 
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where the i sum is restricted to the spatial components q i of the phase space vectors x. 

We next describe the connection between Eq.(1.3) and periodic orbits in the fundamental cell. As the full 
M — >■ M reduction is complicated by the nonabelian nature of G, we first introduce the main ideas in the abelian 
M — >■ M context, and discuss the problems associated with the full reduction in Sect. 3. 



Reduction from M to M 

The general idea of our approach consists in writing 

{e fi<it ~ x) ) M = f ^ dxdy e^*-^ Prob t (x -> y) 

JxeM, y£M 

Here, \M\ = J M dx is the volume of the elementary cell M. Note that there is a unique lattice translation h 
such that y = y — h, with y G M. Therefore, and this is the main point, translational invariance together with 
(1.1) can be used to reduce this average to the elementary cell: 

(e /^-)) M = 1 / dxdy efKW-^y - ^(x)) . (2.1) 

\ M I Jx,y£M 

In this way the global ft flow averages can be computed by following the flow ft restricted to the elementary 
cell M. As is well known [R] , the t — >■ oo limit of such averages can be recovered by means of transfer operators. 
The Equation (2.1) suggests that we study the operator C* whose kernel is given by 

C\y,x) = e?<*'-h(y - z t ) , (2.2) 

where x t = ft(x) G M, but x,x t ,y G M. It is straightforward to check that this operator has the semigroup 
property, J M dz C t2 (y, z)C fl (z,x) = C t2+tl (y,x) . The quantity of interest (1.3) is given by the leading 
eigenvalue of A = e*®^. In particular, for (3 = 0, the operator (2.2) is the Perron-Frobenius operator, 
with the leading eigenvalue A = 1 (the probability conservation). 
To evaluate the spectrum of £, consider 



JM 



(2.3) 



- 1 - 



Here n t (x) is the discrete lattice translation defined in (1.1). For discrete time and hyperbolic dynamics we have 

" £, = £ E | det(1 _j. wr (*•*) 

p:<Tp r =t , x Gp 1 v / / I 

where the sum is over periodic points of all prime cycles p whose period a p divides t, and J p (x) = D^T v (x) . 
Note that the sum over cycle points of p can be replaced by a factor <r p , as det(l — J p ) = det(l — J p (x)) and 
h p = n (Tii (x) are independent of x. For the Jacobian J p this follows by the chain rule, and for the travelled 
distance h p this follows by continuing the path periodically in M. For the discrete time case we finally obtain 

( °° cr p r r/3-hp \ 

-gV |det(i-j;)| j ' (2 - 5) 

where the product runs over the set V of prime cycles. 

Generalization to continuous time [Bo, CE1] amounts to the replacement z CTj> — >■ e~ s ° v , where a p is now 
the (not necessarily integer) period of the prime cycle p: 

zv,.)= nexp -E; | det(1 _ J;) | ) • (*•«) 

p£V \ r=l 1 V p/ I / 

The associated Ruelle ^ function is then (see for ex. ref. [AAC] for details) 

8-h p -s<T p \ 

UC((3,s) = I] 11- , A , , (2.7) 




with A p = Y[ e A the product of the expanding eigenvalues of J p . (This formula can also be obtained by 
using the suspension formula in Appendix C.3 of [R], with the discrete dynamical system leading from collision 
to collision with the central disk, and the ceiling function being the time between these hits.) 

Our first result is therefore: The function Q(/3) of Eq. (1.3) is the largest solution of the equation 
Z((3,Q((3)) = (or equivalently, of 1/(((3,Q((3)) = 0). 

The above infinite products can be rearranged as expansions with improved convergence properties [AAC] . 
To present the result, we define t = e /3 " np ~ S(7p /\A p \, and expand the ( function (2.7) as a formal power series, 

ik i - t p) = i - £' ( 2 - § ) 

PEV Pi,— ,Pk 

where 

ts x = (-l) k t t ■■■t , 

{piviPfej- v ' yi ?2 Pk ' 

and the sum is over all distinct non-repeating combinations of prime cycles. Following the derivation of Eq.(35) 

and Eq.(80) in [AAC] we get, for example, 
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_ E'(-l) k (n pi + --- + n pk )J\K--- A P k \ 
p=o Z'(-mr pi + --- + r pk )/\A pi ---A pk \ 



(2.9) 



with sums as in (2.8). Two derivatives yield our second result: 
The diffusion constant (1.6) is given by 

_ 1 E'(-I) fc (% 1 + --- + %J 2 /|A P1 ---A p j 
2 E'(-l) fc (r pi + --- + r p J/|A pi -.-A p J " 

Note that the global trajectory is in general not periodic, n p ^ 0; nevertheless, the reduction to the elementary 
cell enables us to compute relevant quantities in the usual way, in terms of periodic orbits. 
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Consequences of lattice symmetry 



The lattice symmetry of the Lorentz billiard has important consequences on the properties of the function Q{(3) 
which are best illustrated by introducing its analytic continuation to (3 = ik. The function F(k) = Q(ik) is 
the rate associated with the incoherent scattering function (expzA; • (x t — x)) M considered in light or neutron 
scattering experiments in liquids, in particular, by Van Hove [BY, VH]. The vector k is interpreted as the 
wavenumber of the hydrodynamic modes of diffusion which we also find in the Lorentz gas. The function F(k) 
turns out to be a dispersion relation since F(k) = —Dk 2 + 0(k i ) in an isotropic diffusive system. The isotropy 
of a liquid implies that the dispersion relation only depends on the amplitude \k \ of the wavenumber. 

On the other hand, the lattice symmetry of the Lorentz gas imposes special restrictions on the properties 
of the dispersion relation F(k) and on the values taken by the wavenumber. The present classical problem of 
diffusion is similar to the quantum motion of a particle in a periodic potential. Hence, the wavenumber takes 
its values in the so-called Brillouin zone [BSW, H]. A mode of diffusion is associated with each value of the 
wavenumber k so that the direction of k is privileged in the system. As a consequence, the symmetry of the 
lattice is reduced by the choice of k. This symmetry reduction is formalized by the concept of little group 
associated with the wavenumber which is the subgroup of the lattice point group leaving invariant the vector 
k. For most values of k inside the Brillouin zone, the little group is trivial because it only contains the identity. 
However, the little group is larger when the wavenumber belongs to special symmetry lines or symmetry points 
in the Brillouin zone. In particular, the little group coincides with the full point group when k = 0. 

These results have the following consequences for the factorization of the zeta function Z(k, s). The little 
group associated with k has a single irreducible representation for most values of k where the zeta function does 
not factorize. Only, at those special values of k where the corresponding little group contains several irreducible 
representations, there exists a factorization of the zeta function of the form [CE2] 



where a runs over the irreducible representations of the little group of k. For instance, in the triangular Lorentz 
gas, the Brillouin zone has several lines where the little group contains the identity together with a reflection. 
We can then simplify the calculation of the diffusion coefficient by following the results of Sect. 2 applied to the 
factor of (3.1) corresponding to the identity representation. Nevertheless, the simplification obtained by taking 
advantage of this reflection symmetry is modest since, in this case, the zeta function splits into only two factors. 



Fig. 2. Unfolding of a trajectory of the fundamental domain under the action of the group G. Note that each fundamental domain orbit 
corresponds to 12 distinct global orbits. 



Z(k,s) = 11 Z a (k,s) 



(3.1) 
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The preceding considerations concern the consequences of the lattice symmetry for the factorization of 
the zeta function. There is a different problem which is to express the zeta function in terms of the prime 
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periodic orbits of the fundamental domain M of the lattice (see Fig. 2) rather than those of the elementary 
(Wigner-Seitz) cell M. This problem is a priori independent of the factorization discussed above and presents 
the following difficulty. The stumbling block appears to be the non-commutativity of translations and rotations. 
More precisely, in contrast to (2.4), the global distance 4> raf (x) — 5, x £ p, depends on the starting cycle point 
if p is only a segment of the global cycle p. An example is the diamond- shaped cycle of Fig. 3; depending on 
whether one starts at x 1 or £ 2 , the global distance covered in time <r~ is either the short or the long diagonal. 
This difficulty can be handled in the following way. 




Fig. 3. A fundamental domain 2-cycle p which covers in one return time only 1/2 of the corresponding global cycle p. f(x) is the 
collision-to-collision mapping induced by the flow. To make the figure more readable, the disks have radii smaller than those needed for 
the "no infinite free flight" condition. 

We can introduce the flow </r on the fundamental domain M using the full point group of the lattice. 
Moreover, the density p(x) on which the Perron-Frobenius operator C* acts can always be decomposed using 
the projectors onto the irreducible representations of the full point group. For an arbitrary wavenumber the 
Perron-Frobenius operator will mix the different components of the density p(x), which we can express by 

C\y,x) = R{x;k,t)8[y - 4>*(x)] , (3.2) 

where R(a;; k,t) is a matrix ruling the dynamics of the different components of the density. With Eq.(3.2), the 
Perron-Frobenius operator is reduced to the flow in the fundamental domain M. The zeta function can then be 
written as a product over the prime periodic orbits p of the fundamental domain 
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°° 1 - exp(— sa~r) 

^-tr[R-( — 



Z(k,s) = | | exp^ - X; ~ r tr[R#)1 \ , (3-3) 

where Rp(&) is the matrix R(a;; A;, t) associated with the prime periodic orbitp. It is only when the little group of 
the wavenumber k has nontrivial irreducible representations that the matrices Ilp(k) r split into block-diagonal 
submatrices which can be assigned to each irreducible representation so that the zeta function factorizes as 
explained with Eq.(3.2). The same discussion can be developed with (3 = ik. 

We end this section with the remark that the signature of the lattice symmetry appears in the behavior 
of the function Q(f3 x ,f3 y ) away from f3 = as shown elsewhere [G] and, this, even in triangular or square 
lattices where diffusion is isotropic. Accordingly, the full function contains more information on the 

lattice symmetry than the diffusion matrix of the second derivatives of Q(f3). 



Conclusions 

Compared to the literature [CE2, Ro], the new feature of the problem at hand is use of vector- valued functions 
in Eq. (1.2). The arbitrary vector f3 is only a device for generating moments - the moments themselves are 
invariant under discrete symmetries, but it can be interpreted in terms of the wavenumber of the hydrodynamic 
modes of diffusion as discussed in Sect. 3. 

We have thus obtained a description of global diffusive properties of an infinite periodic dynamical system, 
such as the Lorentz gas, in terms of periodic orbits restricted to the elementary cell. These formulas have been 
tested extensively in refs. [CGS, BEC] on the Lorentz gas, and in ref. [A] on 1 -dimensional mappings. Related 
trace formulas have been independently introduced and tested numerically in ref. [V]. The formalism has been 
generalized to evaluation of power spectra of chaotic time series in ref. [CFP]. However, a derivation of the 
corresponding formulas for dynamics restricted to the fundamental domain should need further developments 
as sketched in Sect. 3. 

In practice, the periodic orbit evaluations of the diffusion constant converge poorly compared with averages 
over scalar quantities such as the Lyapunov exponents. These difficulties are due to several reasons: (1) the 
diffusion coefficient is not a mean but a variance which is always more difficult to evaluate than mean quantities 
like Lyapunov exponents; (2) there is presently no simple formula for the diffusion coefficient in terms of 
the periodic orbits of the fundamental domain; (3) systems like the Lorentz gas do not have simple symbolic 
dynamics and, furthermore, the flow may have discontinuities where trajectories become tangent to the disks. 
These reasons affect the numerical implementation of the formulas proposed here. 
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